New determination of the mass of the rj meson at COSY-ANKE 
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A value for the mass of the r\ meson has been determined at the COSY-ANKE facility through 
the measurement of a set of deuteron laboratory beam momenta and associated 3 He center-of-mass 
momenta in the dp — > 3 He X reaction. The 77 was then identified by the missing-mass peak and the 
production threshold determined. The individual beam momenta were fixed with a relative precision 
of 3 x 10 _J for values around 3 GeV/c by using a polarized deuteron beam and inducing an artificial 
depolarizing spin resonance, which occurs at a well-defined frequency. The final-state momenta 
in the two-body dp — > 3 He?7 reaction were investigated in detail by studying the size of the 3 He 
momentum ellipse with the forward detection system of the ANKE spectrometer. Final alignment 
of the spectrometer for this high precision experiment was achieved through a comprehensive study 
of the 3 He final-state momenta as a function of the center-of-mass angles, taking advantage of 
the full geometrical acceptance. The value obtained for the mass, m v = (547.873 ± 0.005 a tat ± 
0.027 sys t) MeV/c 2 , is consistent and competitive with other recent measurements, in which the 
meson was detected through its decay products. 

PACS numbers: 13.75.-n, 14.40.Be 



I. INTRODUCTION 

The precise value of the mass of the 77 meson has been 
the subject of intense debate for several years. The sit- 
uation seems to have been finally resolved with the pub- 
lication of four experiments that obtained consistent re- 
sults to high accuracy [H-Q. In all these new experi- 
ments the meson was cleanly identified through one or 
more of its decay modes, 7r 7r°7r° [l], 77 [2], or both ||, 
or these plus tt + tt^tt° and n + Tr~j Taking only decay 
experiments into account, the Particle Data Group now 
quote their "best" estimate of the mass as being = 
(547.853 ±0.024) MeV/c 2 ^. Experiments in which the 
77 meson was identified through a missing-mass peak in 
a hadronic production reaction have all reported a lower 
value for the mass, typically by about 0.5 MeV/c 2 . This 
was the case for the ir~p — > nrj reaction, where the beam 
momentum was determined macroscopically to high pre- 
cision using the floating wire technique [6]. In the two 
experiments where the dp (pd) — > 3 Hc 77 reaction was 
used [3, HI, the beam momentum was measured by study- 
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ing other two-body reactions with known final masses. 
One possible concern might be whether the background 
under the 77 peak could have been slightly distorted by 
a strong coupling of, for example, r\ 3 He ^ 7T7t 3 He. Al- 
ternatively, the beam momenta could have been poorly 
determined, though this was done using different tech- 
niques for the three experiments [6|-|8|, or the spectrom- 
eters aligned insufficiently well. 

The situation can only be clarified through the per- 
formance of a much more precise missing-mass exper- 
iment 20]. It is the purpose of the present paper to 
provide results of such an investigation of the two-body 
dp — > 3 He 77 reaction near threshold that is comparable in 
accuracy with those that studied the 77 decay products [3- 
4] . For this to be successfully achieved, it is necessary (i) 
to establish the beam momentum, (ii) to identify well the 
77 meson from the dp — > 3 He X missing-mass peak, and 
(iii) to establish the increase of the final-state momentum 
with excess energy Q = ^/s — (m3 Ho — m v )c 2 above the 
3 He77 reaction threshold with high accuracy. Here y/s 
is the total center-of-mass (cm.) energy. A data set of 
beam momenta and associated final-state momenta near 
threshold then permits the determination of the produc- 
tion threshold and hence the 77 meson mass. 

In order to obtain a clean identification of the 77 me- 
son from a missing-mass peak in a dp 3 He X reac- 
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tion, it is important that rj production be very strong in 
the near-threshold region. Here we are particularly for- 
tunate in that the dp —> 3 Hery total cross section rises 
within the first 1 MeV above threshold to its plateau 
value of rs 400 nb and then remains nearly constant up 
to Q = 100 MeV [H[l3]. This allows one to collect sim- 
ilar statistics over a wide excess energy range, even very 
close to threshold, without expending excessive measur- 
ing time. Furthermore, it has been demonstrated that 
the ANKE facility is well suited for measuring the total 
and differential cross sections of the dp — > 3 He?y reaction. 
The spectrometer has full geometrical acceptance for ex- 
cess energies Q < 15 MeV. The multipion background 
under the rj peak varies smoothly with Q, making a ro- 
bust subtraction of this background possible [tj [l(| . 

The extrapolation of data to identify the production 
threshold requires a precise measurement of the increase 
of the 3 He cm. momentum as a function of the beam 
momentum. Due to resolution or smearing effects, which 
are always present in any real detector, the reconstructed 
momenta can be shifted compared to the true ones. To 
understand and compensate for such effects, a careful 
study of the entire ANKE spectrometer is needed. Af- 
ter calibrating the spectrometer by measuring a variety 
of other nuclear reactions, the requisite precision was 
achieved by demanding that the magnitude of the true 
3 He momentum from the dp — > 3 He rj reaction should be 
identical in all directions in the cm. frame. This check 
was only possible because of the 100% angular acceptance 
of ANKE for the reaction of interest. 

We have previously described in full the measurement 
of the beam momenta via the induced spin depolariza- 
tion technique [ll|. The background subtraction that 
allows the extraction of the rj meson signal from the 
dp — » 3 HcAT reaction is essentially identical to that of 
our earlier work Q. Therefore, most of the emphasis 
here will be on describing the fine calibrations of the 
spectrometer required to get the necessary precision in 
the final 3 He momenta. After outlining the method for 
determining the rj mass in Sec. HH the beam momentum 
determination is briefly summarized in Sec lIIII The stan- 
dard calibration of the ANKE forward detector system 
is given in Sec. IIVI Section [V] discusses how to exploit 
the full geometrical acceptance of a two-body reaction, 
such as dp — > 3 Hc rj, to verify and improve the alignment 
of the spectrometer. Results and estimated uncertain- 
ties are presented and it is shown how resolution effects 
can influence the measurement of the final-state momenta 
and hence the missing mass. 

The results of the determinations of the beam and 
final-state momenta are brought together in Sec. IVII so 
that a reliable extrapolation to the rj threshold can be 
made. In this way, a value for the mass of the rj meson 
was obtained with small systematic and negligible sta- 
tistical errors. The comparison with the results of other 
experiments is made in our conclusions of Sec. I VIII The 
good agreement with other precision measurements was 
only achieved through the exploitation of the large angu- 



lar acceptance of ANKE for this 77-production reaction, 
a feature that was not available in earlier missing-mass 
experiments. 



II. METHOD FOR THE DETERMINATION OF 
THE rj MESON MASS 

By studying the two-body dp — > 3 He?7 reaction in a 
fixed-target experiment, the rj mass can be determined 
by measuring only the momenta of the deuteron beam 
and the recoiling 3 He. In a standard missing-mass exper- 
iment, the rj mass is extracted by measuring the relevant 
kinematic variables at a single fixed energy. This requires 
both a precise calibration of the detector, for a correct 
determination of the 3 He momentum, and an accurate 
measurement of the deuteron beam momentum. 

A much more effective way to measure the rj mass relies 
on the determination of the production threshold by in- 
vestigating the change of the final-state momentum as a 
function of the beam momentum. This requires the kine- 
matics to be measured at several different energies close 
to threshold. For the dp — > 3 He 77 reaction, the final-state 
momentum pf in the cm. frame 
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is a very sensitive function of the rj mass and the total 
energy y/s. The latter is completely fixed by the masses 
of the initial particles and the laboratory momentum, pd, 
of the deuteron beam; 
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The final-state momentum depends only on the 
beam momentum, the 77 mass and other well-measured 
masses [l2j]. If one can fix the production threshold, 
Pf(s) ~ 0, the rj mass can then be determined from 
knowledge of pd- The precision is enhanced because in 
this region dm^/dpd ~ 0.24/c. 

An obvious advantage of the threshold determination 
method is that it does not require a perfect spectrome- 
ter calibration. This is illustrated in Fig. []] using Monte 
Carlo simulated data at twelve different excess energies 
in the range of Q = 1 — 11 MeV and assuming some fixed 
value of the rj mass. Although the analysis was performed 
using Eq. (jTJ) (solid line), to first order p 2 j depends lin- 
early on pd near threshold. 

To reduce the sensitivity of the method to systematic 
errors in the spectrometer calibration, an additional scal- 
ing factor S is introduced which multiplies the right-hand 
side of Eq. ([1]) . Such a factor will occur in a real exper- 
iment, e.g., through minor inaccuracies in the determi- 
nation of the interaction vertex relative to the detection 
system. If one then used data at a single excess energy, 
this would lead to an error in the determination of the 
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FIG. 1: (color online) Identification of the dp — > 3 He?7 pro- 
duction threshold by studying the relationship between the 
square of the final-state momentum p/ and the deuteron beam 
momentum pd- The twelve Monte Carlo points are compared 
with the shape expected from the kinematics (solid line). Al- 
though this shape depends slightly on m v , it is essentially 
a straight line over the range shown. The stability of the 
threshold extrapolation is illustrated by the dashed red and 
blue lines, where the measurements of pj are scaled by arbi- 
trary constants S = 0.9 and 1.1, respectively. 



missing mass. However, by measuring over a range of 
excess energies, such a scaling factor would not affect 
the value of pd corresponding to the production thresh- 
old that is used in the mass determination. This effect 
is illustrated by the dashed lines in Fig. [TJ where 10% 
changes in the scaling factor were considered. 



III. BEAM MOMENTUM DETERMINATION 

A method to determine the beam momentum to high 
precision in a storage ring was developed at the electron- 
positron machine VEPP-2M at Novosibirsk [13] . The 
technique uses the spin dynamics of a polarized beam and 
relies on the fact that depolarizing resonances occur at 
well-defined frequencies, which, apart from the gyromag- 
netic anomaly, depend only upon a particle's speed. By 
depolarizing a polarized beam with an artificial spin res- 
onance induced by a horizontal radio frequency magnetic 
field of a solenoid, the value of the beam momentum can 
be measured to high accuracy. The frequency of such an 
artificial resonance f r is fixed in terms of the machine fre- 
quency fo by the relativistic 7 factor of the particle. Since 
frequencies can be routinely measured with a relative pre- 
cision rs 10~ 5 , an accurate value of 7 and hence momen- 
tum can be deduced. When this method was used for 
the first time at COSY with a vector-polarized deuteron 
beam, the precision achieved was more than an order of 
magnitude better than that of conventional methods [ll[ . 



In the COSY-ANKE experiment, twelve closely-spaced 
beam momenta above the dp — > 3 He r\ threshold were di- 
vided alternately into two so-called supercycles that in- 
volved up to seven different machine settings. In addition 
to six momenta above threshold, one was included below 
threshold in both supercycles to provide the background 
description. Each supercycle covered an excess energy 
range from 1 MeV to ~ 10 MeV. The beam momenta 
were measured both before and after five days of data- 
taking in order to study systematic effects. 

In both supercycles there were collective shifts of the 
spin resonance frequencies of ~ 12—17 Hz, which proba- 
bly originated from changes in the orbit length in COSY 
of about 3 mm. The set of beam momenta from the 
first supercycle were found to decrease slightly over the 
data-taking time whereas those of the second supercycle 
showed the inverse behavior. Average values of the mea- 
surements, before and after data-taking, were calculated 
for each of the twelve beam momenta and these were used 
in the threshold fits. The maximum uncertainty for these 
mean values was conservatively evaluated as ±164 keV/c. 
It is reasonable to assume that the orbits changed linearly 
in time to give a uniform probability distribution of true 
momenta over this interval. The systematic uncertainty 
of the averaged beam momenta is therefore estimated to 
be Apd iSyst = 95 keV/c (rms). Evidence in favor of the 
approach adopted here is found by comparing the results 
obtained with the two supercycles, which are discussed 
in Sec. [VTJ 

The statistical uncertainty in one of the twelve beam 
momenta in the range 3.1-3.2 GeV/c is dominated by 
that of the revolution frequency fo. The measured revo- 
lution frequency deviated randomly by up to ±6 Hz from 
measurement to measurement. The rms uncertainty in 
this frequency, 3.5 Hz, corresponds to a statistical uncer- 
tainty in the beam momentum of Ap<j iSt at = 29 keV/c. 

In total, the twelve beam momenta were measured with 
an overall accuracy of Apd/pd = 3 x 10~ 5 . This is suf- 
ficient to satisfy the needs of a competitive 77 mass mea- 
surement. 



IV. THE ANKE SETUP 



A. The ANKE magnetic spectrometer 



The experiment was performed at the Cooler- 



Synchrotron of the Forschungszentrum Julich 14] using 
the ANKE magnetic spectrometer [n| that is located at 
an internal target station of the storage ring. ANKE 
includes three dipole magnets (Dl, D2, D3) that de- 
flect the COSY beam through a chicane in the ring (see 
Fig. This allows particles produced at small angles, 
which would normally escape down the beam pipe, to be 
detected and this is of particular importance for near- 
threshold reactions. Dl deflects the circulating beam 
by an angle off its straight path onto a cluster-jet tar- 
get [lj| , the spectrometer dipole magnet D2 separates the 
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FIG. 2: (color online) The ANKE setup used for the determination of the 77 mass. The COSY deuteron beam hits the hydrogen 
cluster-jet target and the charged particles produced at small angles are separated by the D2 spectrometer magnet to be 
detected by the ANKE Fd-system. 



produced particles from the beam for momentum anal- 
ysis, and D3 leads the unscattered particles back onto 
the nominal orbit. Although ANKE is equipped with a 
variety of detectors, only the Forward system (Fd) was 
used for the i] mass measurement. This consists of one 
drift chamber, two multiwire proportional chambers, and 
three layers of scintillation hodoscopes [171, [l8| . The lo- 
cation of the various elements of the experimental setup 
is illustrated in Fig. [2] 

The design requires the D2 magnet and the forward 
system to be placed jointly on a moveable platform. The 
deflection angle of the beam, a, depends on the position 
of the platform. The angle (0° < a < 10.6°), the mag- 
netic field strength of D2 (< 1.57 T), and the beam mo- 
mentum cannot be chosen independently of each other. 
By adjusting these three parameters it is possible to in- 
crease the geometrical acceptance. The deflection angle, 
which may differ slightly from the nominal one, was de- 
termined to be 5.8°. 

The tracks of fast charged particles detected in the 
ANKE Fd-system can be traced back through the pre- 
cisely known magnetic field to the interaction point and 
this leads to a momentum reconstruction for registered 
particles. 



B. The standard Fd-system calibration 

In addition to the three parameters mentioned, i.e., the 
magnetic field, deflection angle, and interaction point, 
the drift and wire chamber positions also have to be de- 
termined with high accuracy, and this requires a precise 
calibration of the detection system. In the standard pro- 
cedure [l9|, the positions of the drift and wire chambers 
on the movable platform are first aligned by using data 
taken at the beginning of the beam time with a deflection 
angle of 0° and no magnetic field in D2. The ejectiles then 
move on straight tracks starting from the nominal inter- 
action point in the overlap region of the COSY beam and 



the cluster-jet target. From an analysis of these tracks, 
the position of the Fd-detector relative to the D2 mag- 
net and the target can be reconstructed and compared 
to direct measurements. 

After making this first alignment, the global positions 
of the drift and wire chambers are defined by that of the 
moveable Fd-system platform. Although the positions of 
the platform and interaction vertex are already known 
by direct measurement, a much more precise determina- 
tion of their values is possible by investigating a series of 
reference reactions at all the energies used in the 77-mass 
determination. These are: 

1. Small angle dp — ¥ dp elastic scattering, with the 
fast forward deuteron being detected, 

2. Large angle dp — > dp elastic scattering, with both 
final-state particles being detected, 

3. dp — > ppn charge-exchange, with two fast protons 
being detected, and 

4. dp — > 3 He7r°, with the 3 He nucleus being detected. 

Deuteron-proton elastic scattering in the backward 
hemisphere allows one to verify energy-momentum con- 
servation in the reconstructed four-momenta. For the 
other reactions, the minimization of the deviation of the 
missing mass from the expected value was used to fix 
the positions of the interaction vertex and the Fd-system 
platform. With the parameters extracted at thirteen 
different beam energies, the missing masses were recon- 
structed to an accuracy of « 3 MeV/c 2 for all four chan- 
nels studied. Though good, it is manifestly insufficient 
for a competitive determination of the 77 mass and the 
more refined technique discussed in the following section 
is needed. 
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V. THE 3 He t) FINAL-STATE MOMENTUM 
ANALYSIS 

A. Fine calibration using the kinematics of the 
dp — >• 3 He tj reaction 

The calibration method described in Sec. lIVBl i.e., the 
study of kinematic variables in measured reactions with 
known masses, is standard for magnetic spectrometers 
and was used in such a way in other rj missing-mass ex- 
periments 7, 8]. Unlike these experiments, the ANKE fa- 
cility has full geometrical acceptance for the 3 He from the 
dp — > 3 He ?y reaction up to an excess energy Q ss 15 MeV. 
By taking advantage of this feature and studying the de- 
pendence of the final-state momentum on the 3 He cm. 
angles, the standard calibration could be significantly im- 
proved. 

For a two-body reaction at a fixed center-of-mass en- 
ergy, the final-state momenta pf are distributed on a 
sphere in the cm. frame with a constant radius given 
by Eq. ([1]). We take the z-direction to lie along that of 
the beam, y to be defined by the upward normal to the 
COSY ring, and x = yx z. The relation of Eq. (JTJ can be 
visualized in a simplified way by plotting the magnitude 

of the transverse momentum, p± = yj p\ + p y , versus the 

longitudinal momentum p z for reconstructed events, as 
shown in Fig. [3J For a better visualization of the an- 
gular distribution, each event is weighted with a factor 
l/p_i_. The expected kinematic loci for dp — > 3 He?7 and 
dp —> 3 He ir° are shown by solid lines. It is clear that the 
numbers of events vary along the circles but the method 
used here only depends upon the position of a kinematic 
curve and not on its population. In addition to single me- 
son production, there is a large accumulation of events 
near the forward direction for pj m 350 MeV/c and these 
correspond to two-pion production in the dp — > 3 He(7r7r)° 
reaction. 

The principle of the refined spectrometer calibration 
is the requirement that the momentum sphere should be 
completely symmetric in p x , p y , and p z . It is therefore 
necessary to study the reconstructed momentum p f care- 
fully as a function of the polar and azimuthal angles $ 
and <p. This requires a clean separation of the 3 He 77 sig- 
nal from the background and this is the first step in the 
analysis. 

In the event selection, 3 He particles were identified 
and the raw background, consisting mainly of deuterons 
and protons from dp elastic scattering and protons from 
deuteron breakup, was suppressed by cuts on the energy 
loss and time of flight of charged particles in the Fd detec- 
tor. Figure. UJa) shows the energy loss times the square of 
the particle velocity versus the laboratory momentum. A 
projection of these events onto the ordinate gives a clear 
peak with a low background. In our analysis we used a 
very loose cut at the 6cr level but, as will be shown in 
Sec. I VII changing this to 3er would have only a tiny effect 
on the value obtained for m„. 




p z [GeV/c] 

FIG. 3: (color online) The magnitude of the reconstructed 
transverse cm. momentum p± in the dp — > 3 HeX reaction 
plotted against the longitudinal cm. component p z at an 
excess energy Q — 6.3 MeV with respect to the 77 threshold. 
For better visualization of the angular distribution, each event 
is weighted with a factor l/p±. The small and large circles 
correspond to the kinematic loci for the dp — » 3 He 77 and 
dp — > 3 He7r° reactions, respectively. ANKE covers the full 
solid angle for r\ production near threshold whereas, for pions, 
only the forward 3 He are detected. 



Figure. 0|b) shows a typical example for one scintil- 
lation counter combination of the length of the recon- 
structed flight path of the 3 He selected from the events 
in Fig. 0Ja). This has been calculated from the TOF 
information between the first and the last scintillation 
wall in the forward detector and the reconstructed par- 
ticle momentum, assuming the mass and charge of the 
3 He. A clear peak is again evident with only a moderate 
background. Reducing the 3cr cut used here to 2a gives 
an even smaller change in the value of m n . 

The remaining background, shown in Figs. [3] and [5l 
originates mainly from some residual deuteron breakup 
and multipion production in the dp — > 3 HeA" reaction, 
where X = (tttt) or X — (tttttt) . At the lowest ex- 
cess energy, the signal/background ratio is around 11 but 
this decreases with increasing excess energy to w 1.8 at 
Q = 10.4 MeV. This background was subtracted using 
the data taken below the 77 threshold at an excess energy 
of Q ss —5 MeV. These data were kinematically trans- 
formed to positive Q in order to compare them with re- 
sults obtained above threshold. The details of this tech- 
nique are described for missing-mass spectra in Ref. [§], 
but the method is equally applicable to final-state mo- 
mentum spectra. Due to the very high statistics of the 
current experiment, the distribution inpf could be inves- 
tigated for twenty bins each in 7? and <p. This is illustrated 
in Fig. [SJ where examples of the pf spectra summed over 
(f> are shown for six cost? bins for the energy closest to 
threshold, Q = 1.1 MeV. A similar picture is found for 
the 4> dependence after summing over 9. 

Mean values of the 3 He momentum pf and the peak 
widths for different 7? and <f) bins were extracted from 
the background-subtracted dp — > 3 He77 distributions by 
making Gaussian fits. A variation of the width of 4- 
12 MeV/c (rms) was found, as well as a displacement of 
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FIG. 5: Center-of-mass distributions of the 3 He momentum 
from the dp — » 3 HeX reaction for six typical polar angle 
bins at the lowest excess energy measured, Q — 1.1 MeV. 
The experimental data summed over <j> are shown by black 
lines and the background estimated from subthreshold data 
by gray lines. The resulting background-subtracted dp — > 
3 He r\ signal is shaded gray. 



FIG. 4: (color online) (a) Two-dimensional distribution of 
energy loss times the square of the particle velocity versus its 
laboratory momentum. A projection of this onto the t/-axis 
gives a peak with little background. Data within a ±6<r cut 
were retained for further analysis, (b) A typical flight length 
distribution for events selected from panel (a). Data within a 
±3cr range were used in the subsequent analysis. 



the mean value, both of which depended upon the polar 
and azimuthal angles. This striking effect results from 
the different resolutions of the ANKE Fd-system in p x , 
p v and p z , which are discussed in the following section. 

B. Influence of the momentum resolution on the 
reconstructed Anal-state momentum 

The influence of resolution on the angular dependence 
of the reconstructed 3 He momentum and the missing- 
mass distributions for the dp — > 3 He rj reaction is illus- 
trated by the two-dimensional (p±,p z ) projection shown 
in Fig. Ufa) . In the ideal case of a measurement with per- 
fect resolution, the final-state momenta are distributed 
on a sphere of constant radius pf, as indicated by the 
black line. Both the missing mass and final-state momen- 
tum are then independent of cost? and <fi, as illustrated 
by the black line in Figs. [5£b)-6(e). 

In a real experiment, the reconstructed momenta in 



the laboratory frame are smeared by the finite resolu- 
tion associated with the detector setup and reconstruc- 
tion algorithms. If only p z were smeared, with a Gaussian 
width of say a z — 30 MeV/c in the laboratory frame, an 
event on the momentum sphere, indicated by the black 
arrow in Fig.[5]Ja), could be shifted along the red arrows. 
It is important to note that for p™ c < pj ue events are 
shifted toward lower | cos i9| whereas the reverse is true 
for p™ c > pj uc . This effect leads to a pf distribution 
that is a function of cos$ (Fig. [S^b), red circles), i.e., 
the momentum sphere is stretched for large longitudinal 
momenta and compressed for high transverse momenta. 
For simple kinematic reasons, the missing mass shows the 
inverse behavior [see Fig. HJd)]. Since the smearing was 
here assumed to be independent of p x andp y , and hence 
4>, this is reflected in the constancy of the reconstructed 
momentum in Fig[B]Jc). Nevertheless, its value is higher 
than the true one, pj uc . 

When only the transverse momenta are smeared, 
but with different Gaussian widths, e.g., (<j x , a y , a z ) — 
(10, 20, 0) MeV/c, the reconstructed momentum has the 
opposite dependence on cosd. As shown by the blue 
crosses in Fig. [Hb), Pf decreases for cos$ w ±1 and in- 
creases for cos & ~ 0. Different resolutions in p x and p y 
also introduce a dependence on 0, and this leads to os- 
cillations in both pf c and the missing mass in the plots 
of Figs. [3c) and 6(e). The amplitude and phase of these 
oscillations depend on the ratio cj x /a y . 
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FIG. 6: (color online) Influence of resolution on the determi- 
nation of the final-state momentum. The ideal pf sphere of 
panel (a) (black) is changed along the horizontal (red) arrow 
by finite resolution in the longitudinal (z) direction, whereas 
resolution effects in the transverse direction are indicated by 
the vertical (blue) arrow. Panels (b) - (e) show the possible 
distortions at Q — 1.1 MeV, evaluated in Monte Carlo simula- 
tion. The mean values of both the final-state momentum and 
the missing-mass distributions for individual cos i? and (f> bins 
are shown without (black line) and with momentum smearing 
in the ^-direction (red circles) and transversely (blue crosses). 



In reality, all three momentum components are recon- 
structed with finite and generally different resolutions 
and the effects described above will be superimposed, 
though they will be dominated by the component with 
the worst resolution. The determination of the rj mass 
has to take these kinematic resolution effects into account 
because, without so doing, the value extracted for m v 
would depend on the production angle. For the current 
ANKE experiment, differences in m n of up to 0.5 MeV/c 2 
are found between cosi? = ±1 and cosi? = 0. Further- 
more, the average of the final-state momentum over all 
cos •& and <j> is shifted to a higher value than the true one 
and the missing mass shifted to a lower value. The angu- 
lar distribution of the dp — » 3 He r\ reaction could modify 
slightly the effects of the resolution but, since this vari- 
ation is linear in cosi? over all the Q-range studied 
even this is of little consequence for the determination of 
rrirj, given the symmetry in cost? shown in Fig. (S^b). 



Precise determinations of the resolutions in (p x ,p y ,p z ) 
are absolutely essential in order to correct the measured 
kinematic variables. How this is done, by demanding an 
isotropic pf momentum sphere, will be discussed next in 
some detail in Sec. IV CI 



C. Correction to the final-state momenta 

Mean values of the measured final-state momenta for 
background-subtracted dp — > 3 He 77 distributions are 
shown in Fig. [7] for twenty individual cos $ and 4> bins 
at an excess energy Q — 1.1 MeV, before and after the 
improvement of the calibration. The results for the stan- 
dard calibration, presented in the upper panels, show 
that the momentum sphere is neither centered nor sym- 
metric. The momentum sphere is shifted to higher p z , 
i.e., on average pf is higher for 3 He produced in the for- 
ward direction than in the backward. The oscillations 
in the <j> spectrum are also far from being symmetric, 
and this is particularly evident at <f> ps ±90°, where the 
p y momentum component dominates. This asymmetric 
pattern is rather similar at all twelve energies and this 
stresses the need to improve the calibration for the de- 
termination of the correct momenta. 
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FIG. 7: (color online) The mean values of the final-state mo- 
mentum distributions are shown for individual cosi? and <j> 
bins for the standard (top) and improved (bottom) calibra- 
tion at Q = 1.1 MeV (red circles). The results of Monte 
Carlo simulations are shown without (black horizontal line) 
and with momentum smearing (black points). The compari- 
son of the data with the simulation leads to a determination 
of the momentum resolutions in the three directions. 

The p y component depends sensitively on the relative 
y-position between the two wire chambers in the forward 
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detector and, by varying their positions by about 0.3 mm, 
the momentum sphere could be centered in p y . Changes 
in the magnetic field strength of 0.0015 T at 1.4 T, i.e., 
changes of the order of 0.1%, allow one to make the neces- 
sary adjustment in the p z component. The p x component 
can be fine-tuned by adjusting the deflection angle a by 
about 0.4% at 5.8°. By making changes such that the 
mean values of the final-state momenta are distributed 
on a centered and perfectly symmetric sphere in cos *d 
and <j), the detector alignment can be significantly im- 
proved, as shown in the lower part of Fig. [7J This proce- 
dure was carried out using the data at all twelve energies 
simultaneously. The magnitudes of these changes are so 
small that they have no impact on the values of the miss- 
ing masses of the different reactions used in the standard 
calibration. 

The improved spectra, shown in the lower half of Fig. [7] 
for one of the twelve energies, allow one to study the mo- 
mentum smearing in the three directions. The values 
of <r x and cr y were determined from the amplitude and 
phase of the oscillation in (j>, as explained in Sec. IV Bl 
whereas that of a z was extracted by making a second 
order fit to the data and simulations for pf as a func- 
tion of cos-d. An additional constraint is that the val- 
ues of (o Xl a y ,cr z ) must reproduce the width of the 3 He?/ 
final-state momentum signal when integrated over all $ 
and 4>. The resolution parameters were determined from 
the spectra with uncertainties of (Aa x , Aa yi Aa z ) = 
(0.2, 0.2, 0.1) MeV/c. 

The individual momentum spreads were determined 
separately for each of the twelve energies. In contrast to 
a Zl which is constant to within 1 MeV/c, a x decreases 
and a v increases with excess energy. This behavior is 
reasonable because the cone angle of the 3 He ejectilcs 
from the dp — » 3 He r\ reaction increases with excess en- 
ergy. Close to threshold the hit positions in the cham- 
bers are located in a small area near the center of the 
wire chambers, whereas at higher Q the hits are more 
widely distributed. With the resolution parameters thus 
determined, data similar to those shown in Fig. [7] are de- 
scribed well by Monte Carlo simulations at all twelve en- 
ergies. Figure. [5] shows the distribution of the final-state 
momentum summed over all cos i? and <fi at an excess en- 
ergy of Q « 8.6 MeV. It is clear that the measured mo- 
mentum distribution (filled area) is generally very well 
reproduced by the Monte Carlo simulated data (black 
crosses) . 

Using data from all twelve energies above threshold, 
a mean momentum resolution of ANKE was calculated 
in the narrow 3 He laboratory momentum range of 2.63- 
2.68 GeV/c. The values of the momentum spreads in 
the laboratory frame were found to be (<r x , a y , a z ) — 
(2.8, 7.9, 16.4) MeV/c. As expected for a fixed target 
experiment, the resolution in p z is by far the poorest. 
Furthermore, because of the particular construction of 
the wire chambers, the p x resolution is better than that 
for p y . However, in the determination of the final-state 
momentum corrections, the individual resolution param- 
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FIG. 8: (color online) Final state momentum distribution for 
background-subtracted dp — > 3 He r\ data (filled red area) and 
simulation (black crosses) for all events at an excess energy of 
Q ~ 8.6 MeV. The vertical line indicates the true final-state 
momentum assumed in the simulation. 



eters determined at each beam momentum were used. 

Owing to the resolution effects shown in the lower half 
of Fig. [?l the average of the final-state momentum over 
all cos and <f> is shifted to a higher value than the 
true one (black horizontal line). By comparing the av- 
erages resulting from the Monte Carlo simulations with 
and without momentum smearing, correction parameters 
were calculated for all twelve energies and these are pre- 
sented in Fig. [S] The correction is about 2.22 MeV/c for 
the lowest momentum and decreases steadily with pf. 
The error bars shown are dominated by the uncertain- 
ties in the resolution parameters {Aa x , Aa y , Aa z ) and 
range from 0.08 MeV/c at the lowest beam momentum 
to 0.04 MeV/c at higher energies. It should be noted 
that the dependence of the correction parameters on the 
value assumed for the rj mass is negligible. 

If the resolution factors <Ji are largely independent of 
the beam momentum, the correction should vary like 
r~> 1/pf. This behavior arises because the deviation de- 
pends on the ratio of the ANKE momentum resolution to 
the size of the momentum sphere. Confirmation of such a 
dependence is offered by the curve, which is a 1/pj fit to 
the data. Despite the good x 2 /NDF ps 0.9, the individ- 
ual Monte Carlo estimates were used in the subsequent 
analysis. 

Compensation for the effects of the momentum res- 
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FIG. 9: (color online) Deviation of the measured final-state 
momentum in the dp —> 3 He?7 reaction from the true one 
due to resolution effects, evaluated in Monte Carlo simula- 
tion. The twelve measured final-state momenta in the near- 
threshold region, Q = f-II MeV, have to be corrected by 
0.7-2.2 MeV/c to compensate for such effects. The curve is a 
1/pf fit to the points. 



olution is essential for an accurate determination of 
the production threshold. Without this correction, the 
value obtained for the rj mass would be lower by about 
150 keV/c 2 . 



VI. THE MASS OF THE rj MESON 

In order to obtain a robust value for the mass of the 
•q meson, it is necessary to extrapolate the experimental 
data with high precision in order to determine the value 
of the deuteron beam momentum at threshold. For this 
purpose we show the complete set of final-state and beam 
momenta (pf,pd) with statistical uncertainties in Table[TJ 
Also shown are the corrections that were included in the 
values quoted for p/. 

Due to the resolution effects, the twelve pf distribu- 
tions are found to be slightly asymmetric over the whole 
angular range, both for the experimental data as well as 
the Monte Carlo simulated events, as shown for a typi- 
cal energy in Fig. [5] The values of the pf means were 
determined for both data and simulations within ±3cr 
limits. The good statistics of « 1.3 x 10 5 3 Her] events 
for each energy meant that the uncorrected value for pf 
could be extracted with an uncertainty of rs 23 keV/c. 
In total w 1.5 x 10 6 3 He?y events were collected in the 
experiment. After applying the resolution correction to 
the measured pj, the overall uncertainty is increased, as 
indicated in Table [TJ 

Although the description of the reconstructed final- 
state momentum distribution in Fig. [8] by the Monte 
Carlo simulation is very good, it is not perfect, espe- 
cially in the high momentum tail. Such discrepancies 
could arise from slight imperfections in the spectrometer 



calibration, 3 He scattering in the wire chambers, or limi- 
tations in the background subtraction approach. In order 
to quantify their influence on the value extracted for the r\ 
mass, the interval used to determine the means of the pf 
was varied between ±2cr to ±4er, where a represents the 
peak width. Such a variation leads to a collective shift 
in the extracted final-state momenta of approximately 
0.16 MeV/c. Since this effect corresponds to an over- 
all shift in the final state momenta, it is not included in 
the numbers quoted in Table Q] but must be considered 
in the final rj mass determination, where it introduces a 
systematic uncertainty of 12 keV/c 2 . 



TABLE I: Values of the laboratory beam momenta pd, cor- 
rected final-state cm. momenta pf, and the pf correction 
parameters measured at twelve different excess energies; the 
statistical uncertainties are noted in brackets. The approx- 
imate values of Q quoted here are merely used to label the 
twelve settings. 



Excess 


Beam 


Final-state 


Pf correction 


energy Q momentum momentum pf 


parameter 


MeV 


MeV/c 


MeV/c 


MeV/c 


1.1 


3146.41(3) 


32.46(8) 


2.22(8) 


1.4 


3147.35(3) 


35.56(7) 


1.94(6) 


1.6 


3148.45(3) 


39.00(6) 


1.69(6) 


2.1 


3150.42(3) 


44.09(6) 


1.57(5) 


2.6 


3152.45(3) 


49.25(5) 


1.35(5) 


3.1 


3154.49(3) 


53.66(5) 


1.28(4) 


4.1 


3158.71(3) 


61.70(5) 


1.07(5) 


5.1 


3162.78(3) 


68.77(4) 


0.95(4) 


6.3 


3168.05(3) 


76.92(4) 


0.88(3) 


7.3 


3172.15(3) 


82.64(5) 


0.84(4) 


8.6 


3177.51(3) 


89.81(4) 


0.76(4) 


10.4 


3184.87(3) 


98.64(4) 


0.74(4) 



The extrapolation of the data to threshold is illustrated 
in Fig. [10] for both pf and pj versus pd- Whereas, to 
first order, pj depends linearly on pd, the analysis con- 
siders the full dependence Pf — pf(m v ,S,Pct), as given 
by Eqs. (UJ and $2$. Only the r\ mass, chosen as a free 
parameter, defines the production threshold. The scal- 
ing factor S, discussed in Sec. |TTJ allows for a possible 
systematic energy dependence of pf. This would repre- 
sent yet a further fine tuning of the description of the 
measurement process but it is crucial to note that its 
introduction does not affect the value obtained for the 
threshold momentum and hence m n . 

The overall fit to the data in Fig. [TU] has a x 2 / NDF = 
1.28 and the best value of the mass quoted in Table ITU 
is m v = (547.873 ± 0.005) MeV/c 2 , where the error is 
primarily statistical. The corresponding deuteron mo- 
mentum at threshold is p d = (3141.688 ± 0.021) MeV/c. 
A linear fit of pj versus pd would give a poorer reduced 
X and a mass that was 10 keV/c 2 higher. 

The scaling factor S — 1.008±0.001 is well determined 
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beam momentum p rf [GeV/c] 

FIG. 10: (color online) Corrected values of the final-state cm. 
momentum p/ (black crosses) and its square (red stars) plot- 
ted against the deuteron laboratory momentum pd- The error 
bars are too small to be shown on the figure. The extrapola- 
tion to threshold is carried out on the basis of Eq. where a 
scaling factor S has been introduced. The lower panel shows 
the deviations of the experimental data from the fitted curve 
in pj. The errors shown here do not include the overall sys- 
tematic uncertainty in pf associated with the description of 
the profile in Fig. 



TABLE II: Values of the 77-mass and scaling factor evaluated 
separately for the two supercycles and for the complete data 
set. The errors do not include the systematic uncertainties in 
the determination of the beam momentum. 

Supercycle Scaling factor S m v 

1 1.008 ± 0.001 (547.870 ± 0.007) MeV/c 2 

2 1.008 ± 0.001 (547.877 ± 0.007) MeV/c 2 
1+2 1.008 ± 0.001 (547.873 ± 0.005) MeV/c 2 



and differs very slightly from unity, which means that the 
twelve momentum spheres are about 0.8% bigger than ex- 
pected. One consequence of this is that the missing mass 
is not constant and one would get a slightly different 
value at each of the twelve energies studied. If one took 
S = 1.0, the x 2 /NDF would jump to 24.7 and this would 
result in a shift of 64 keV/c 2 in the 77 mass. Such a sys- 
tematic error is avoided in the threshold-determination 
method described in Sec. [TTJ 

By far the dominant systematic errors arise from the 
determinations of the absolute value of the beam momen- 
tum and the pf correction parameters. As can be seen 
from Table UTTl all other sources, such as effects from the 
time stability of the data, further contributions from the 
fine calibration, the event selection, the background sub- 



traction for the pf distributions, as well as contributions 
of the t] mass assumed in Monte Carlo simulations, are 
negligible in comparison. 

TABLE III: Systematic uncertainties in the determination of 
m v . The small "experimental settings" contribution includes 
effects from the magnetic field, the deflection angle, and the 
(vertical) wire chamber positions, all of which are coupled. 
The PDG value of m, [5|] was used in the simulations but, 
if our result were used, it would only result in a 2 keV/c 2 
change. The effects of putting stricter cuts on AE x f3 2 and 
the flight length are also shown. 



Source 


Variation 


Am, 
keV/c 2 


Absolute beam momentum 


95 keV/c 


23 


Experimental settings 






2 


m v assumed in simulations 


20 keV/c 2 


< 2 


AE x p 2 cut 


6(7 - 


•>■ 2<T 


5 


Flight length cut 


3a - 


■¥ 2a 


1 


Pf correction parameters 


4<r - 


* 2a 


12 


Total systematic uncertainty 






27 



This uncertainty in the beam momentum translates 
into one in the mass of 

Am, = -. mpPd - A Pd = 23 keV/c 2 , (3) 

(mane + rn v )E d 

and hence, taken together with all other systematic un- 
certainties, to a final value of 

m v = (547.873 ± 0.005 sta t ± 0.027 syst ) MeV/c 2 . (4) 

To investigate further some of the systematic effects, 
the results were extrapolated separately for the data ob- 
tained in the two supercycles and the individual values of 
the r\ mass and the scaling factor S are given in Table QTJ 
There is only a tiny difference between the two separately 
determined r\ mass values of 7 keV/c 2 . This agreement 
supports the validity of taking the mean values of the 
beam momenta determined at the beginning and end of 
the measurements and the subsequent correct handling 
of the corresponding systematic uncertainties. 

VII. CONCLUSIONS 

We have measured the mass of the 77 meson in a 
missing-mass experiment by identifying precisely the pro- 
duction threshold in the dp —> 3 Her] reaction. As is seen 
in Fig. [Til the value obtained is consistent with all the re- 
cent measurements where the meson decay products were 
studied 0-01 • The precision achieved is similar to these 
works and the deviation from the PDG best value [|[ is 
only 20 keV/c 2 , which is less than our systematic error. 

This success was based upon a precise determination 
of the beam momentum using the spin-resonance tech- 
nique [11| , a clear identification of the rj signal [9j , and a 
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FIG. 11: Results of the different r\ mass experiments. Where 
two error bars are shown, the heavy line indicates the statisti- 
cal uncertainty and the faint ones the systematic. The earlier 
missing- mass experiments, marked Rutherford Lab. 74 [f|, 
SATURNE 92 0], and COSY-GEM 05 [|, all obtained low 
values of m n compared to the experiments where the meson 
was identified through its decay products, viz. NA48 02 [H, 
KLOE 07 0, CLEO 07 H, and MAMI-CB 12 [|. Our re- 
sult, COSY-ANKE 12, is completely consistent with these 
more refined experiments. 



rather than yields to determine the threshold momentum 
and hence the meson mass. However, it is important to 
realize that the anomalous behavior of the production 
cross section, where the cross section jumps so rapidly 
with excess energy, leads to the desirable high count rates 
near threshold. 

Our result differs by about 0.5 MeV/c 2 from earlier 
missing-mass evaluations [6|-[8j and so the hypothesis of 
a distortion of the background under the r\ peak must 
be discarded. Unlike experiments with external targets, 
the energy loss in a windowless cluster jet is negligible, 
but the current experiment had other advantages over 
these measurements. In particular, if we had access only 
to data taken in the forward direction, this would not 
have allowed the fine tuning of the spectrometer and a 
somewhat different value would have been found for m,. 

Finally, twelve energies above threshold were inves- 
tigated and this allowed a reliable extrapolation to 
be made to find the production threshold. This is 
intrinsically subject to far fewer systematic uncertainties 
than a measurement at a single energy. It is therefore 
clear that, with care, a missing-mass approach can be 
competitive with experiments in which meson decays 
are measured. 



systematic study of the measurement of the 3 He7? final- 
state momentum in the ANKE spectrometer. The latter 
was made possible only through the complete geometric 
acceptance of ANKE for the dp — > 3 He r\ reaction close to 
threshold. This allowed us to require that the cm. mo- 
mentum in the final state should be identical in all direc- 
tions. This is a powerful technique that might be useful 
for other two-body reactions. Unlike the MAMI method- 
ology the experiment relied purely upon kinematics 
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